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^ , ABSTRACT 

■ We use the galaxy st ellar mass and hal o merger tree information from the semi-analytic model 

galaxy catalogue of iFont et alj (120081) to examine the accretion of galaxies into a large sam- 

Q . pie of groups and clusters, covering a wide range in halo mass (10 12 9 to 10 153 h~ x M Q ), 

' and selected from each of four redshift epochs (z=0, 0.5, 1.0 and 1.5). We find that clusters 

ryj \ at all examined redshifts have accreted a significant fraction of their final galaxy populations 

through galaxy groups. A 10 14 5 hr x M Q mass cluster at z=0 has, on average, accreted ~ 40% 
of its galaxies (M ste u a r > 10 9 M Q ) from halos with masses greater than 10 13 hr 1 M Q . 
Further, the galaxies which are accreted through groups are more massive, on average, than 
galaxies accreted through smaller halos or from the field population. We find that at a given 
' epoch, the fraction of galaxies accreted from isolated environments is independent of the final 

lO i cluster or group mass. In contrast, we find that observing a cluster of the same halo mass at 

■ each redshift epoch implies different accretion rates of isolated galaxies, from 5-6 % per Gyr 

at z=0 to 15% per Gyr at z=1.5. We find that combining the existence of a Butcher Oemler 

qq ■ effect at z=0.5 and the observations that galaxies within groups display significant environ- 

<T^> ' mental effects with galaxy accretion histories justifies striking conclusions. Namely, that the 

Q\ , dominant environmental process must begin to occur in halos of 10 12 - 10 13 hr 1 M Q , and 

■ act over timescales of > 2 Gyrs. This argues in favor of a mechanism like "strangulation", 

in which the hot halo of a galaxy is stripped upon infalling into a more massive halo . This 
simple model predicts that by z=1.5 galaxy groups and clusters will display little to no envi- 

■ ronmental effects. This conclusion may limit the effectiveness of red sequence cluster finding 

methods at high redshift. 
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1 INTRODUCTION 

In recent years, an extraordinary confluence of independent mea- 
surements of the cosmological parameters has resulted in the con- 
cordance model of the Universe (ACDM), in which the mass den- 
sity is dominated by cold dark matter. In this model, the ini- 
tial distribution of density perturbations has the greatest power 
on small scales, which causes low mass dark matter haloes to 
form first at high redshift. Larger haloes form later through the 
merging, or accretion, of smaller halos. Eventually, this 'hier- 
archical structure formation' leads to the formation of galaxy 
groups and clusters, which become more common with time. The 
mass assembly of dark matter halos has been extensively stud- 
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ied analytically (Press & Schech teil[l974l;lBond et alJll99ll;lBowe 
1 199 it iLacev & Coldl 19931; ISheth & Tormenlboojlvan den Bo~ 
2002; Benson et al. 2005) and through numerical simulations 
(Davis et al j 1 19851; Li et alj 120081) . Consistent with these studies, 
Berrie ret al J ( 2009L hereafter B09) used n-body simulations to 
show that the mass assembly of clusters is dominated by the most 
massive accretion events; in effect, the merging of groups with clus- 
ters. However, by associating dark matter subhalos with galaxies, 
they show that the galaxy assembly of clusters is dominated by 
lower mass halos, or the infalling of isolated galaxies. This dis- 
tinction could be of great importance since there are a variety of 
physical processes that depend on the mass of the host dark matter 
halo and which could affect the properties of a galaxy, such as ram 
pressure stripping, strangulation and galaxy harassment. 

Indeed, detailed observations of dense environments, galaxy 
groups, and clusters in the local universe have shown that the galax- 
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ies which inhabit these environments have properties substantially 
different from galaxies in low density or field environments. In 
particular, galaxy gro ups and clusters have lower average galaxy 
star formation rates dLewis et alj|2002l; [Gomez et alj|2003 | ), lowe r 
fractions of disk galaxies (Dressier et al. 1997; McGee et al. 2008), 
and higher red fractions dBalogh et al.l2004l ; IWeinmann et al . 2006) 
than field galaxies. Despite this wealth of observational data, there 
is no consensus on the dominant physical mechanism responsible 
for these differences, mainly because large populations of "tran- 
sition" objects have avoided detection. In particular, there is no 
large excess in the fraction of galaxies betwe en the red sequence 
and the blue cloud i n dense environments ( Balog h"et alj 120041 ; 
Weinmann et al ] |2006l) . While there are specific examples of tran- 
sitioning spiral galaxies which are in the process o f having their 
HI gas stripped due to ram pressure in local clusters dKennev et al.l 
2004; IVollmer et alj |2004), the X-ray temperatures and pressures, 
as well as the infalling velocity of the galaxies, required for such a 
transformation mechanism are probably too high to be effective in 
low mass groups. 

Strangulation, the process in which the more loosely bound 
hot halo of a galaxy is stripped by the group or cluster halo, 
leaving a reduced amount of gas available for future star for- 
mation (Baloghetal. 2000), is an attrac tive candidate because 
it is still effective in low mass groups ^McCarthy etail l2008bl : 
iKawata & Mulch aev 2008). However, it is not clear if such a gen- 
tle mechanism can account for t he dramatic effect seen in clus- 
ters. |z^iudoff^Muicha^ ( 1 19981) have proposed that the extreme 
properties of galaxy clusters may result from the "pre-processing" 
of galaxies in group environments before accretion into the cluster. 
This is supported by observations of reduced star for mation rates 
in the outskirts of cluster s, well past the virial radius (Balo gh et al.l 
Il999l iLewis et al.ll2003) . However, B09 have claimed that "pre- 
processing" is not a large effect. They find only ~ 12 % of galax- 
ies are accreted in to the final cluster environment as members of 
groups with five or more galaxies. While the B09 clusters are rela- 
tively low mass, their work shows the importance of distinguishing 
the accretion of galaxies from that of dark matter mass. 

A complementary approach to trying to isolate "transition 
galaxies" is to study the properties of galaxies in groups and clus- 
ters as a function of redshift. As first shown by Butch er & Oemlerl 
ill 9781) and confirmed by many others (eg. [Laverv & Henrvl[l986l ; 
ICouch & Sharpies! 1987l : [Ellingson et al.l200ll) . the fraction of blue 
galaxies in clusters increases with redshift, the so called Butcher- 
Oemler effect. Despite this, the fraction of star forming galaxies 
in groups and clusters is still lower than the coeval field fraction 
at lea st to z=l dWilman et all2005l ; lGerke et all2007r . lBalogh et al.l 
2009). The need to explain the Butcher-Oemler effect, as well as the 
local properties of galaxy clusters provides important constraints 
for the nature of the transformation mechanism. Essentially, if the 
transformation mechanism only occurs in very massive clusters, 
then the fraction of blue galaxies is simply the fraction of galaxies 
which have fallen into the cluster within the time scale of transfor- 
mation. 

The time scale for transformation of galaxy properties to oc- 
cur is a significant uncertainty in attempting to link the growth 
of structure to the Butcher-Oemler effect. Previous attempts us- 
ing cluster assembly histories adopted relatively short time scales 
of ~ 1 Gyr and, while complicated by uncertain cosmological pa- 
rameters, showed that a direct infall mo del alone did not produce 
enough evolution in the blue fraction dBowed[l99ll ; iKauffmannl 
1995). iKodama & Bower (200 ll) combined the evolving star for- 
mation properties of field galaxies with a cluster infall model to 



successfully reproduce d the scatter in the red sequence of low red- 
shift clusters. Similarly, lEllingsonet"ai]d200ll) found that the radial 
distribution of early type galaxies in galaxy clusters at two redshift 
epochs could best be explained if the galaxy infall into clusters de- 
creased by a factor of ~ 3 between z > 0.8 and z ~ 0.5. 

In this paper, we examine the galaxy assembly properties of 
groups and clusters over a wide mass range and at four redshift 
epochs. We investigate the mass of halos through which groups and 
clusters gain their galaxies and the extent to which preprocessing 
in the group environment is important at four redshift epochs. By 
making simple assumptions, we investigate the predictions for the 
fraction of galaxies in groups and clusters which are "environmen- 
tally affected" for a range of relevant timescales and the halo mass 
thresholds which those effects begin. Using these models we try to 
gain insight into the dominant physical processes necessary to re- 
produce observations of group and cluster galaxies, as well make 
predictions for future observations. In Sj2] we present the details of 
our simulated clusters and some of their properties and present our 
results in Sj3] We discuss these results and conclude in Sj4] In this 
paper, we assume a cosmology with Q m = 0.25, S1a = 0.75, 08=0.9 
and Hq = 100ft km s^ 1 Mpc -1 unless mentioned otherwise. 



2 SIMULATIONS 

To interpret observations of galaxy properties as a function of envi- 
ronment, we need to know the accretion history of those galaxies; 
as shown by B09 this can be subtly different from the total mass 
assembly history. But galaxy formation has proven to be a diffi- 
cult problem, and it is not clear, given that the dark matter halo 
mass function has a very different shape from the galaxy luminos- 
ity function, if the approach of B09 of simply associating subhalos 
with galaxies includes all of the relevant physics. At the least, this 
approach does not allow for the robust identification of the stel- 
lar masses of galaxies. Unfortunately, an obvious alternative — the 
direct simulation of the baryonic processes of galaxies — is diffi- 
cult on the scale of the cosmological volumes needed to study large 
samples of groups and clusters. 

Semi-analytic galaxy formation models provide a good tool to 
encapsulate the essential p hysical processes of gas cooling, star for- 



mation and feedback (e.g. [ White & Frenk 1991; Kauffmann et al 
1993 |; ISomerville & Primackll999l ; ICroton et all2006l ; lBower et al 



2006). Dark matter simulations, on which modern semi-analytic 
models are based, are now large enough to allow the study of the 
growth of the groups and clusters over a wide r ange of redshifts . 
We make use of one such semi-analytic model bv lFont et al.l d2008l . 
hereafter F08), which is a recen t modification to th e Durham semi- 
analytic model (GALFORM) o flBower et all d2006h . The basic pre- 
scriptions for gas coo ling and star for mation in the GALFORM 
model was laid out bv lColeetall d2000l) . an d subsequently modi - 
fied for mod ern cosmolog i cal pa rameters bv lBenson et alj f2003l) . 
The model o f lBoweretal.ld2006l) introduced a method for parame- 
terizing the effect of AGN feedback on the gas in massive galaxies 
to correc t for the "ove r coolin g" problem. 

The iBoweTe t al. (2006|) model, as in essentially all previous 
semi-analytic models, implements a relatively simple treatment of 
environmental effects, in which the hot gas reservoirs of galaxies 
are removed upon becoming a satellite galaxy. Many authors have 
since shown that this approach produces an unphysically high frac- 
tion of red galaxies in groups and clusters ( Weinman n et alj|2006l ; 
iBaldrv et alj|2006l ; iGilbank & Baloghll2008l) . The F08 model im- 
plements a more realistic "strangulation" model in which the hot 
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gas halo of galaxies falling in to more massive halos a re removed 
according to a prescription of iMcCarthv et af (2008b). However, 
a careful examination of cluster and group data with this model at 
a range of redshifts reveal that there are important discrepancies. 
In particular, the model overpo pulates the green v alley between 
the blue cloud and red sequence Balogh et al] d2009l McGee et al., 
in prep.). We emphasize that despite this difficulty in reproducing 
galaxy colours, the stellar masses of galaxies in the F08 and Bower 
et al. models are much better understood. In particular, the Bower 
et al. model reproduces the observed evolution of the stellar mass 
function out to at least z=5. 

In this paper, our analysis relies primarily on the GALFORM 
prediction of the stellar mass function of galaxies in different envi- 
ronments. This is insensitive to the problem noted above, as the star 
formation rate of galaxies declines rapidly with redshift, so the bulk 
of a galaxy's stellar mass is already in place before it ever becomes 
a satellite. Thus, the details of the strangulation procedure adopted 
in GALFORM are unimportant for our analysis and, indeed, all our 
conclusions are independent of the choice of either the Bower et al. 
model or the F08. model. 
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2.1 Cluster and group sample 

The F08 model, from which our simulated galaxy clusters and 
groups are drawn, is based on merger trees derived from the 
dark matter Millennium simulation ( S pringel et al .l2005t) , a ACDM 
cosmological box with 50 0//i Mpc sides . The Millennium sim- 
ulation uses GADGET2 JSpringelll20bl) . a TREE-PM N-body 
code, and an initial power s pectrum calculated using CMBFAST 
(Seljak & Zaldarriaga 1996). Th e merger trees are generated as de- 
scribed in Hellv et all J2003h and lHarker et alj J20061) . and are com- 



plete down to halos which host ~ 10 8 h^ 1 Mq galaxies. In this pa- 
per, we are principally concerned with selecting samples of galax- 
ies which are observationally accessible, and thus specify a single 
fixed stellar mass cut of M > 10 9 h~ x Mq, much higher than the 
completeness limit. 

We analyze all the groups and clusters in the F08 model more 
massive than M = 10 12 ' 9 hT 1 Mq at four redshift epochs (z=0, 
0.5, 1 and 1.5). The key properties of the cluster samples are shown 
in Table Q] In particular, we show the number of clusters, and the 
average number of galaxies with stellar masses above M = 10 9 
h~ Mq at the epoch of observation, in each of the mass bins which 
will be used in the remainder of the paper. 

In Figure [T] we present the cumulative distribution of galax- 
ies which reside within the virial radius of host halos of a given 
mass. We plot this for four stellar mass ranges at z=0. In the F08 
model, ~ 50% of z = L» galaxies are in host halos with masses 
above 10 12 ' 5 ft" 1 M ra . This compares very well with observational 
results: iBerlind et all f20o3) found that ~ 56 % of M r < -20.5 
galaxies in the SDSS are linked to groups containing at least one 
other member, a result that is co mpletely c onsistent with indepen- 
dent analysis using the 2dFGRS lEke et all We also see that 25% 
of L„ galaxies are in relatively large groups or clusters with halo 
masses above 10 13 h~ x Mq. This is much larger than the ~ 10 
% claimed by B09, likely a result of the way they assign galaxies 
to subhalos, as discussed further in <j3 .61 In particular, B09 assign 
a galaxy to a subhalo if the subhalo mass is > 10 1 h~ Mq 
when it is accreted into a more massive host. However, the mass in 
a subhalo begins to be tidally stripped significantly before reaching 
the virial radius of a more m assive host, even with out significantly 
disturbing the galaxy within (Natarajan et al. 20071). 



Figure 1. The cumulative distribution of the host halo mass of galaxies at 
z=0. The distribution is shown for four ranges in the galaxy's stellar mass 
at z = 0, shown in the upper right corner in units of h~ 1 Mq. 



3 RESULTS 

We now look in detail at how the cluster galaxies end up in the 
clusters, and what insights this might give into the processes which 
might affect those galaxies. 



3.1 Cluster and group accretion history 

Galaxies which have been in massive halos prior to joining the final 
environment may have been environmentally pre-processed. Thus, 
we begin by examining the host halo masses of galaxies just prior 
to their accretion into the final group or cluster halo. To acheive 
this, we trace the most massive progenitor of every galaxy, back 
through each snapshot in the simulation. We record the halo mass 
of this progenitor in the timestep just before it becomes a member 
of the final cluster, which defines its environment at the time of 
accretion. 

We show the full accretion histories for all the cluster mass 
bins, in each of the four redshift epochs, in the Appendix. Here we 
will examine the most important insights which can be drawn from 
those accretion histories. Figure [2] shows the fraction of galaxies 
in the final cluster which were accreted through haloes at least as 
massive as 10 13 hT 1 Mq (large groups). We show this as a function 
of the final cluster mass for each of four redshift epochs. We first 
consider relatively low-mass clusters, with M ~ 10 14 ' 2 Mq at z=0. 
We find that 32 per cent of galaxies in these clusters were accreted 
through such group-sized halos. This is somewhat higher than the 
24 per cent found by B09; the small difference can be related to 
the difference in the way haloes are populated with galaxies, as we 
discuss in § 13.61 

However, such clusters are fairly poor systems; they are less 
massive than all 16 clusters observ ed extensively by the CNOC1 
collaboration fcarlbergetaflll996l) . and an order of magnitude 
smaller th an the nearby Com a cluster (M 2 oo = 1-88 j^gjj X 10 15 
fr 1 Mq, dKubo et alj[2007T) ). Figured shows that the fraction of 
galaxies which are accreted through group sized halos is strongly 
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Redshift 
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40 


15.0-15.6 


15.14 


1161 
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14.7-15.0 


14.82 


569 
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14.53 
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13.94 
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13.8-14.1 
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13.62 
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Table 1. Properties of the cluster sample derived from lFont et alj (2008). The first column lists the redshift snapshot from which the clusters were selected 
and the second column gives the total number of clusters used for analysis in each bin. Columns 3 and 4 list the cluster halo mass range and median mass of 
clusters in that range. We use these halo mass bins extensively in the rest of the paper. Column 5 lists the average number of galaxies per cluster with stellar 
masses above M = 10 9 h~ 1 Mq at the epoch of observation. 



dependent on the mass of the final halo. This is because massive 
haloes are not surrounded by an average patch of the unive rse, but 
tend to be strongly clustered with other massive halos (eg. iKaisej 
1 19841) . At z=0, we see that 45 % of galaxies accreted into a cluster 
of Coma's mass have been accreted from haloes with M > 10 13 
h^ 1 Mq. This suggests that pre-processing in group environments 
before cluster accretion may be significant. Interestingly, the frac- 
tion of galaxies accreted through massive haloes has only a weak 
dependence on the redshift of observation of the cluster. In other 
words, a Coma-sized cluster at z=0.5 would accrete 40 % of its 
galaxies from M > 10 13 h -1 Mq halos. The galaxy assembly his- 
tories are remarkably similar, with the dominant difference being 
simply that Coma-sized clusters do not exist in the relatively small 
volume of the Millennium simulation at z = 1.0 or 1.5. 

In Figure [3] we show the fraction of stellar mass which is ac- 
creted through halos at least as massive as 10 13 h' 1 Mq. This 
figure is quite similar to Figure [2] However, notice that the frac- 
tion of stellar mass accreted by the most massive clusters through 
groups is larger than the fraction of galaxies accreted through such 
systems. Indeed, the stellar mass accretion history closely matches 
the expected behavior of the dark matter accretion. The extended 



Press Schechter formalism and n-body simulations of dark matter 
roughly agree that ~ 30 % of the mass of a halo is accreted from ha- 
los with masses a tenth the mass of the final halo teond et aljl99ll : 
lBowerll99lMLacev & Coleil993l ; IStewart et alj2008l) . We find this 
same fraction for all our stellar mass accretion histories, while the 
fraction of galaxies accreted is smaller at high cluster mass. This 
implies there are are fundamental differences in how galaxies are 
accreted as a function of their stellar mass. This is illustrated in 
Figure [4] where we show the accretion histories of galaxies which 
end up in aM = 10 15 h~ x Mq cluster at z=0, binned by their 
final stellar mass. There is a large difference in the masses of the 
host halos prior to accretion for low and high mass galaxies. While 
~ 52 % of the most massive galaxies are accumulated from haloes 
with M > 10 13 hT 1 Mq, this is only the case for ~ 45 % of the 
least massive galaxies we consider. This discrepancy is much larger 
if we consider accretion through poorer groups, with M > 10 12 
h^ 1 Mq. The more massive galaxies are more likely to have been 
accumulated from group mass halos, and thus more likely to have 
been pre-processed prior to accretion into a cluster. 

Now that we have seen that the degree of group preprocess- 
ing depends on both the stellar mass of the galaxy and the mass 
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Figure 2. The fraction of cluster galaxies which were accreted into the final 
cluster halo as a member of a halo with M > 10 13 h~ 1 Mq . This is shown 
as a function of the final cluster mass at the epoch of observation, for four 
redshifts. All cluster galaxies have final stellar masses of M > 10 9 Mq. 
The mass range bins were defined in Table 1, and are shown for all bins 
containing more than two clusters. 



Figure 4. The cumulative distribution of accreted cluster galaxies as a func- 
tion of host halo mass at the time of accretion into the final cluster. The dis- 
tribution is shown in three stellar mass bins at z=0, for a final cluster with 

M = lO 15 ^- 1 Mq. 
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Figure 3. As Figure [2] but showing the fraction of accreted stellar mass 
which resides in a M > 10 13 h^ 1 Mq halo at the time of accretion. 



of the final cluster, we would like to examine how this varies as 
a function of redshift. In Figure [5] we show the fraction of cluster 
galaxies which were accreted into the final cluster halo as a mem- 
ber of a halo with M > 10 13 h~ x Mq halo. This is broken up into 
three bins, which represent the redshift at the time of the galaxy's 
accretion into the cluster. From this we see that the degree of pre- 
processing is significantly dependent on the time the galaxies were 
accreted. Galaxies which are accreted recently into the cluster are 




13.5 14 14.5 15 

Mass of Cluster [Log(M halo (fr 1 M s ))] 

Figure 5. The fraction of cluster galaxies which were accreted into the final 
cluster halo as a member of a 10 13 h~ 1 Mq halo or greater. This is shown 
as a function of the final cluster mass z=0 and for three bins in accretion 
redshift. All cluster galaxies have final stellar masses of M > 10 9 Af Q . 



more likely to have been in a group environment than ones accreted 
into the cluster at high redshift. In particular, since 2 = 0.5 the most 
massive clusters today have accreted most of their new galaxies via 
infalling groups. 
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3.2 Cluster and group assembly histories 

We have seen that the accretion history of clusters varies with final 
cluster mass, is a function of the stellar mass of the accreted galaxy 
and is dependent on the redshift of accretion. However, this does 
not address the state of the cluster itself. The importance of pre- 
processing depends not only on the accretion history but also on 
the amount of time the main cluster progenitor itself had the mass 
of a group. 

Therefore, to get a complete picture of the assembly of galaxy 
clusters and groups and the halo masses which are important for the 
properties of their galaxies we present Figure[6] This shows the dis- 
tribution of halo masses in which the most massive progenitors of 
final z=0 cluster galaxies reside, as a function of lookback time and 
for four bins of final z=0 cluster mass. The panels in Figure[6]show 
distinctly different assembly histories for very massive clusters, 
smallish clusters, and groups. In particular, the relative importance 
of the group environment varies tremendously for these three types 
of structures. The most massive cluster never has more than 17 % 
of galaxies in group sized halos (10 13 hT 1 M < M na io < 10 14 
h^ 1 M ) while as many as 44% of the galaxies in a M ~ 10 14 ' 2 
h^ 1 M© cluster today have spent some time within such haloes in 
the past. In fact, for a period of 2 Gyrs, groups are the most com- 
mon environment of the galaxy progenitors; this is because during 
this time the main cluster progenitor itself has the mass of a group. 
Thus, considering only the haloes of galaxies prior to accretion into 
the main cluster may underestimate the role of the group environ- 
ment, as already noted by B09. 

Given the distinctly different assembly histories of these clus- 
ters and massive groups, it is perhaps surprising that observations of 
large samples of galaxy clusters in the local universe show that the 
fraction of red galaxies is app roximately constant in clusters more 
massive than 10 13 ' 8 h' 1 M ^Hansen et ai]|2007h . Therefore, it is 
useful to look for some common trait in the assembly histories of 
clusters which may point to the reason for this uniformity. It is in- 
teresting that the population of 'isolated' galaxies, those in 10 11 
h^ 1 M < Mhaio < 10 12 ft -1 M©, shows a similar distribution in 
the four different panels. At a lookback time of 10 Gyrs, ~ 55 % of 
cluster galaxy progenitors were in this halo mass regime, and that 
percentage has declined at a nearly constant rate of 5 - 6 % per Gyr 
until the current epoch, regardless of the final cluster mass. In other 
words, the distribution of galaxies not in 'isolated' halos is similar 
regardless of final cluster mass, and supports the hypothesis that the 
galaxy transformation mechanism begins to occur as galaxies leave 
their 'isolated' halos. 

Finally, we examine the assembly histories of galaxy clusters 
of a given mass at each redshift epoch. In Figure [7j we show the 
distributions of halo masses for the most massive progenitors of 
10 14 ' 5 ft -1 M cluster galaxies as a function of lookback time at 
all four redshift epochs. While the final cluster mass is the same (at 
each epoch), the higher redshift clusters must assemble their mass 
more quickly and thus their galaxies have not been in massive halos 
for as long. For instance, 5 Gyrs prior to the observation epoch, ~ 
50% of z=0 galaxies were in 10 14 M h' 1 haloes, while none of 
the z=l or z=1.5 cluster galaxies were even in 10 13 ft" 1 M haloes 
yet. Environmental processes have had a much longer timescale 
over which to affect low redshift groups and clusters than higher 
redshift ones. 

There are two interesting points when comparing Figure UJ 
with Figure [6] First, we see that the maximum fraction of galaxies 
in each halo mass bin is the same in clusters of the same final mass 
but seen at different redshifts. For instance, the maximum fraction 



of galaxies which reside in halos of 10 h~ Mg< Mhaio < 10 
h^ 1 Mq at any time is 40 % regardless of the redshift epoch. The 
lookback time at which these maximum fractions occur varies sig- 
nificantly with redshift, but it would appear their path through the 
hierarchy is similar. Essentially, clusters of fixed mass at differ- 
ent redshift epochs have assembly histories which become more 
stretched out at lower redshift. The assembly histories would look 
almost identical if the lookback time was divided by the age of the 
universe at that redshift epoch. This result was hinted at in Figure 
[2] which showed that the fraction of galaxies accreted through mas- 
sive haloes was approximately the same at all redshift epochs for a 
cluster of given mass. 

This leads to the second interesting observation to be made 
from Figure[7j The rate at which galaxies leave their 'isolated' ha- 
los increases significantly with redshift. At z=0, as before, for the 
10 Gyrs prior to observation the fraction of galaxies in halos of 
10 11 h' 1 M < M ha i < 10 12 h' 1 M decreases by about 5-6% 
per Gyr, while 10% (15%) [20%] of galaxies leave their 'isolated' 
halos per Gyr at a constant rate for 5(3.5)[2.5] Gyrs prior to ob- 
servation at z=0.5(l)[1.5]. Therefore, the accretion rate of galax- 
ies from isolated environments into groups and clusters is higher 
at higher redshift. Again, this result is a direct result of the re- 
duced time between the epoch of observation and the beginning 
of the universe. The assembly histories at higher redshift are just 
compressed, leading to a higher accretion rate, even though the to- 
tal number accreted from isolated environments is constant at each 
epoch of observation. The effect this has on the galaxy properties 
of galaxy clusters as a function of redshift will be discussed in the 
following section. 

3.3 Cluster to cluster variation in environmental effects 

We have established the galaxy accretion history and galaxy as- 
sembly history of galaxy clusters at a range of epochs. We would 
now like to assess how these galaxy histories affect the final prop- 
erties of galaxies at each redshift epoch. To this end, we exam- 
ine the fraction of galaxies in each cluster, which have been within 
dense environments long enough to expect that environmental ef- 
fects might be important. By examining the fraction of environ- 
mentally affected galaxies in each cluster, we can quantify both the 
total numbers of affected galaxies, and their variation from cluster 
to cluster. 

In a simple way, we can parametrize the length of time it 
takes for a galaxy to display an environmental effect, Tt runc , after 
falling into a halo with a mass above a characteristic mass thresh- 
old, Mtrunc- Although it is not obvious that there is a single main 
physical mechanism which causes the environmental affects dis- 
played in both groups and clusters, we explore the predictions of 
such a model and discuss the limitations of this approach in the 
following section. 

Given this model we can explore how varying the truncation 
time, Tt r u n c, and the characteristic mass threshold, Mtrunc, alters 
the implied environmental effects on galaxies. In Figure(8] we show 
the predicted average fraction of galaxies in each cluster which are 
subject to environmental effects in our simple model; the distribu- 
tion of this fraction is reflected in the four contour lines marking 
the 10, 33, 67 and 90 percentiles. Here we fix the truncation time, 
Ttrunc, to be 3 Gyr and allow the characteristic mass threshold, 
Mtrunc, to vary from 10 12 h' 1 M to 1O 14 M . In other words, in 
this figure, a galaxy has felt an "environmental effect" if it has been 
within a halo of mass M > Mtrunc for at least 3 Gyrs. In addition, 
we allow a fourth category, in which the expression of an environ- 
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Figure 6. The fraction of galaxies residing in z=0 clusters that are found in halos of mass M at a previous time t. We only consider the most massive progenitor 
of each cluster galaxy. Each curve shows a range of M, indicated by the legend, in units of h~ 1 Mq. Each panel represents clusters of different final masses, 
as indicated in the top-left corner. Note that the curve corresponding to haloes with 10 11 < M/h" 1 Mq< 10 12 is very similar in all cases, indicating that 
the accretion rate of "isolated" galaxies is roughly independent of final cluster mass. 
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Figure 7. As in Figure|6] but for a cluster of mass 10 14 5 (h 1 Mq) observed at four redshift epochs as indicated at the top of each panel. Clusters of a given 
mass at higher redshift must assemble their mass more quickly, and thus the time available for pre-processing through group-sized haloes is decreased. 



mental effect occurs 3 Gyrs after the galaxy has become a satellite 
galaxy in a larger dark matter halo, regardless of its mass. 

This figure has some noteworthy features. First, for massive 
clusters (M > 10 14 5 h' 1 Mq) at z=0, the mean number of en- 
vironmentally affected galaxies in this model is similar (~ 80-85 
%) regardless of Mtrunc- The implication of this for low redshift 
observational studies is that it is difficult to discern the value of 
the characteristic mass threshold by observing systems above that 
mass threshold. This highlights the importance of studies of low 
mass galaxy groups. Observations at low and intermediate redshift 
show that group galaxies with a given stellar mass have proper- 
ties distinct from the average field galaxy; if our simple model of 
environment-driven transformation is correct, this indicates a char- 
acteristic mass threshold of at least this scale (M as 10 12 5 - 10 13 
ft -1 M Q ) dWilman et alJl2005l ; IWeinmann et alj20o3) . 

Indeed, as previously mentioned, low redshift observations 
show that the fraction of red galaxies in clust ers is essentially uni - 
form, for clusters with M > 10 13 ' 8 h' 1 M dHansen et al1l2007h . 
Given this, it is also worth noting that in Figure [8] our model 
also produces a strikingly flat fraction of environmentally affected 
galaxies per cluster as a function of cluster mass. This is a direct re- 
sult of the behavior noted in Figure [6] that the fraction of galaxies 
infalling from isolated halos is independent of halo mass. 



Although it may be difficult to use the average properties 
of massive clusters at a given epoch to discern the characteris- 
tic mass threshold, one possible method would be to observe the 
variation in their properties. The predicted scatter in the fraction 
of environmentally affected galaxies per cluster is quite small (~ 
5 %) for 10 14 ' 5 h' 1 M Q clusters at z=0 when M trun c= 10 12 
tr 1 M Q , but close to 40 % when Af t runc= 10 14 ft" 1 Mq. The 
scatter in, for instance, the fraction of early type galaxies or opti- 
cal lin e emitting galaxies in clusters at z=0 is much smalle r than 
40 % jDressledll980l : IPoggianti et alj|2006t iFinn et alj|2008l) . We 
will examine the scatter in red fractions of galaxies in clusters at 
z=0 in a future paper. Unfortunately, the scatter at z=0 of a model 
where M trU nc= 10 12 h~ x Mq is not that different from a model 
where Af t runc= 10 13 hT 1 Mq. However, notice that the scatter in 
these two models bec omes more significant at z > 0. Intriguingly, 
iDressler et al.l dl997t) showed that, while the morphology-density 
relation was equally strong in all clusters at low redshift, the rela- 
tion was stronger in centrally-concentrated clusters than irregular 
clusters at z ~ 0.50 In effect, this suggests that the scatter in the 
fraction of environmentally affected galaxies of each cluster is sig- 



1 While the Dressier et al. results, and many intermediate redshift results, 
have limiting stellar mass on the order of 10 10 h~ 1 Mq compared with our 



8 McGee et al. 



1.0 

0.8 
0.6 

.E 0.2 
k 0.0 



O 0.8 



O 
_Q 
CO 



Q. 



0.6 
0.4 



C/) 
CO 

ns 0.2 
E 
o 



jS i.o 

.g o.8 

£ 0.6 
00 

=5 0.4 



0.2 
0.0 



O 1.0 
CO 0.8 
£ 0.6 
0.4 
0.2 
0.0 




1.0 -z = 0.5, T = 3Gyr, M = 12- 



0.0 ■_ 



— I 1 1 1 1— 



0, T = 3 Gyr, M = 12-- z = 0, T = 3 Gyr, M = 13 



H 1 1 1 1- 



-z = 1.0, T = 3Gyr, M = 12- 



-i — 1 1 — 1 H H i — 

i = 1.5,T = 3Gyr, IVI = 12 



— I 1 1 1 1 — 



z = 0.5, T = 3 Gyr, M = 13-- z = 0.5, T = 3 Gyr, M = 14-- z = 0.5, T = 3 Gyr, M = sat 



r 



-H 1 1 1 h- 



-H 1 1 1 1- 



z = 1.5, T = 3 Gyr, M = 13-- 1.5, T = 3 Gyr, M = 14 



1 1 1 1 1 

z = 0, T = 3Gyr, M= 14 



-H 1 1 1 1- 



z = 1.0, T = 3 Gyr, M = 13-- z= 1.0, T = 3 Gyr, M = 14-- 1.0, T= 3 Gyr, M = sat ■ 



H 1 1 1 1- 



— i 1 1 1 1 — 

z = 0, T = 3 Gyr, M = sat ■ 



h 1 1 1 1- 



h 1 1 1 1- 



z = 1 .5, T = 3 Gyr, M = sat 



_l I l_ 



13.0 14.0 15.0 13.0 14.0 15.0 13.0 14.0 15.0 13.0 14.0 15.0 

Halo mass of Cluster [Log(M SO | ar /h)] 



Figure 8. The fraction of cluster galaxies with M > 10 9 /i~ x Mq that have resided within a halo of mass M >Mt runc for a time t >Tt rU nc is shown as a 
function of final cluster mass. We interpret this as the fraction of "environmentally-affected" population in our simple model. The panels contain four contour 
lines marking the 10, 33, 67 and 90 percentiles of the distribution in this fraction, while the dashed yellow line represents the average. The truncation time is 
fixed at Tt r unc= 3 Gyr, and each row shows a different assumption for Aftrunc> as indicated. Different rows correspond to clusters at a different redshift, as 
indicated. 



nificant at z ~ 0.5. Although not definitive, this may point to a 
characteristic mass threshold which is somewhat larger than 10 12 
h~ x m q , given that scatter in that model is still quite small at z=0.5 
(~ 13 % at 10 14,5 h' 1 M ). Notice that a model where the envi- 
ronmental effects begin to occur when a galaxy becomes a satellite 
behaves very similarly to a model with AI tTunc = 10 12 h^ 1 Mq . We 
discuss this similarity further in t]3.7l 

Examining the redshift evolution of any of the given models 
shows that they all predict a significant Butcher-Oemler effect. That 
is, they predict that there are fewer environmentally affected galax- 

limit of 10 9 h" 1 Mq, we have verified that the scatter in the cluster red 
fractions is constant with a limiting mass change to 10 10 h _1 Mq. 



ies in clusters with increasing redshift. In particular, by z= 1.5 all 
of the models predict a very small or non-existent fraction of en- 
vironmentally affected galaxies. Indeed, the 10 14 h^ 1 M model 
leads to the prediction that, by z — 1, no galaxies will be environ- 
mentally affected. 

Our choice of r trunc = 3 Gyrs in the models presented above is 
ad hoc, and we would like to quantify how changing the timescale 
effects the predictions. In Figure [9] we explore a model in which 
the characteristic halo mass, Mt run c, is kept fixed at 10 12 h' 1 M , 
and allow Tt runc to vary from 1 Gyr to 4 Gyrs. We show the fraction 
of environmentally affected galaxies for each of the four redshift 
epochs of our clusters. Although Mtrunc is held constant, we note 
that the results and our interpretation are similar for any choice 
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Figure 9. As Figure[8] but where the characteristic halo mass threshold is fixed at Aftrunc= 10 12 h 1 Mq, and the truncation times Tt runc are varied along 
rows of the figure, from 1—4 Gyr as indicated. 



of A/trunc within the range 10 12 -10 13 h' 1 Mq, which seems the 
most likely value given the arguments above. 

Similarly to Figure [8] for each Ttrunc we see a significant 
Butcher-Oemler effect, such that clusters at higher redshift have 
fewer galaxies affected by environmental processes. However, the 
size of the effect even between z = and z — 0.5 is dramatically 
altered by the choice of time scale. With a short timescale of only 1 
Gyr, the fraction of environmentally-affected galaxies evolves little, 
from ~85% at z = 0.5 to ~ 95% today. On the other hand, a long 
timescale of Ttrunc = 4 Gyr results in a much stronger evolution 
over this redshift range, from 50% to 80%. Compare this evolution 
with that observed in the red fraction of cluster galaxies, which in- 
dicate an evolution of ~ 25 % over a similar reds hift range, from 
0.9 at z=0.2 to 0.65 at z=0.5 tellingson etal]|200lh . This seems to 
indicate that a relatively long time scale for the expression of envi- 
ronmental effects (> 2 Gyr) would be required to match this quick 



evolution. A similar timescale is necessary to explain the radial gra- 
dient of passive ga laxies in galaxy clusters jEllingson et~aT]|200ll; 
iBalogh et all200(j|) . 

The predicted scatter from cluster to cluster is also notewor- 
thy. Recall that in Figure [8] we saw that the scatter was sensitive 
to the characteristic halo mass used. In this plot, for the majority 
of the time, the scatter is similar at each redshift regardless of the 
timescale for truncation. This strengthens our previous argument 
that a well-defined measure of the scatter in cluster properties at 
a given redshift could allow one to discern the characteristic halo 
mass for truncation. 

We have provided strong evidence, which we summarize in 
£13.51 that the dominant environmental processes at work in galaxy 
groups and clusters begin to become effective at a halo mass scale 
of 10 12 - 10 13 h -1 Mq, and are active for a timescale of at least a 
few Gyrs. Given these constraints, we see that figure[9]predicts that 
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Figure 10. The fraction of environmentally affected galaxies in clusters at 
all four redshift epochs as a function of galaxy stellar mass. Cluster galaxies 
have host halo masses greater than 10 14 h _1 M0at the epoch of observa- 
tion. This model assumes Mt rmlc =10 12 ft _1 Mq and Tt rU nc= 3 Gyrs as 
in Figure [8] 

by z=1.5 there should be little to no environmental effect on galax- 
ies. Remarkably, this p rediction has some o bservational evidence 
to suggest it is correct. ICooper et alj d2007l) showed, using galax- 
ies selected from the DEEP2 redshift survey, that the red fraction 
only weakly correlates with overdensity at z ~ 1.3. While the com- 
parison to our predictions is complicated because the Cooper et al. 
sample only includes massive galaxies, this is not a trivial agree- 
ment; in fact, assuming the timescale was 1 Gyr, this would lead us 
to predict that 70 % of galaxies at z=1.5 are still environmentally af- 
fected. This fraction would be even higher when we used the same 
limiting stellar mass as Cooper et al., as is discussed in the next sec- 
tion. Additionally, the DEEP2 survey is complicated by their rest 
frame blue magnitude limit which causes them to naturally detect 
fewer and fewer red galaxies at higher redshifts. Further, DEEP2 
does not cover a wide enough area to have massive clusters within 
it, so targeted, stellar mass limited studies of the extreme cluster 
environments are still needed at this redshift to quantify the size of 
the environmental effects. Although, given our results, they may be 
difficult to find using the popular and e fficient red sequence method 
dGladders & Yeell200rj|:lLu et alj2009l) . 

3.4 Stellar mass dependence of environmental effects 

Observations suggest that the fraction of galaxies w hich are pas- 
sive or red, depends grea tly on their own stellar mass dBaldry et al.l 
l200dlHaines et alj|2006l) . It is thought that this is at least partially 
due to secular influences , ie. AGN feedback, wh ich primarily oc- 
cur in massive galaxies dKauffmann et alj|2003l) . We can use our 
simple model for environmental effects to examine the fraction of 
cluster galaxies (those with Mhaio > 10 14 ft -1 Mq) which may 
also be subject to environmental effects. This is presented in Figure 
[l0]for a model which has M trU nc= 10 12 hT 1 M Q and T trunc = 3 
Gyr at all four redshift epochs. 

The fraction of environmentally affected galaxies is a strong 



function of stellar mass in this model, and the gradient becomes 
stronger with increasing redshift. The most massive galaxies have 
resided within group-sized haloes since at least z = 1.5; thus any 
environmental effects would have manifested themselves a long 
time prior to observation, and we expect to see little signature of 
cluster growth in their properties. On the other hand, galaxies with 
lower stellar mass are a better tracer of the recent mass accretion 
history of the cluster, hence we see a strong evolution in the frac- 
tion of environmentally-affected galaxies. 

3.5 Observational constraints 

It is now useful to review the observational constraints on our 
model parameters, Mtmncand T trun c- The halo mass threshold at 
which environmental effects become important must be at least as 
low as 10 13 Mq because there are observatio ns of systems 
at thi s mass with significant environmental effects (Wil manet al.l 
120051 ; IWeinmann et alj|2006l) . For this reason, we investigated a 
model with a low halo mass threshold, 10 12 /i" 1 Mq, in Figure|9] 
Well defined samples of galaxy clusters show a significant Butcher- 
Oemler effect, such that the fracti on of red galaxies decr eases from 
~ 0.9 at z=0.2 to ~ 0.65 at z=0.5 dEllingson et alj200ll) . This evo- 
lution is much quicker than predicted by a model with a T trunc of 2 
Gyrs or less. Thus a model with A/t run cOf ~ 10 12 h' 1 Mq and a 
Ttruncof ~ 3 Gyrs is the most favored model. As suggested previ- 
ously, this leads to the prediction that by z=1.5, little or no environ- 
mental effects are felt by the galaxy population. 

Recall £|3,4| in which we investigated the stellar mass de- 
pendence of the galaxy population using our most favored model. 
We found that while the most massive galaxies are environmen- 
tally affected at all redshifts, the lower mass g alaxies become more 
affected with time. iGilbank & Baloghl {2008) used a compilation 
of the observational literature to show that the ratio of red bright 
galaxies to red faint galaxies steadily increases with redshift, the 
same qualitative behaviour we see in the simple model. 

It is difficult to observationally quantify the extent to which 
massive galaxies are environmentally affected. This is largely be- 
cause the visual colours of galaxies are not very sensitive to low 
levels of star formation. Mid-IR observations are more sensitive to 
low levels of star formation and thus are better at establishing the 
environmental influence of massive galaxies. Observations at z ~ 
0.4 suggest that only 10 % of massive galaxies (> 10 10 M ) 
in groups have IR emission indica tive of activity, whil e the global 
frac tion is much high er (~ 40%, Iwilman et alj|2008l) . Addition- 
ally, |WolfetalJ 12009) find that massive galaxies are uniformly old 
and red in the cluster cores, while having a significant population of 
dusty, star-forming red galaxies in the infall regions. Both of these 
studies suggests that significant environmental effects are felt even 
by massive galaxies, as assumed in our model. 

We emphasize that this observational comparison is qualita- 
tive, yet highly suggestive. In a future paper we investigate the 
quantitative behavior of these models with a direct comparison to 
the best available cluster, group and field data to z~ 1 . 

3.6 Comparison to Previous Work 

In an attempt to explain observations of the f raction of cluste r mem- 
bers with [Oil] emission at z=0 and z=0.6, Pogaianti et al. ( 2006) 
have presented a similar, but more complex model. The observa- 
tions they present (their Figure 4) show that while, at z=0.6, higher 
mass clusters have lower average fractions of [Oil] emitting galax- 
ies, this is largely because of an upper envelope which decreases 
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with increasing cluster velocity dispersion. In contrast, they notice 
that at z=0, the fraction of [Oil] emitting galaxies is constant with 
cluster velocity dispersion above 550 km/s (~ 10 14 ft -1 M©), but 
the scatter is large below that value. 

In effect, to explain the observations, Poggianti et alj J2006l) 
introduces two Mttuiw and two Ttrunc parameters to match the ob- 
served behavior. The first set of iWtiunc and T tIunc are meant to rep- 
resent 'primordially' passive galaxies, and are associated with el- 
liptical galaxies. They claim that galaxies within 3x 10 12 Mq 
groups at z = 2.5 represent these primordially passive galaxies. 
The second set of parameters are associated with quenched galax- 
ies or SO galaxies, and are set to have M tr unc=10 14 Mq and 
Ttrunc= 3 Gyrs. However, observations of galaxy groups with 
masses less than 10 14 h -1 MqsIiow a significant populati on of SO 
galaxies and passive spiral galaxies JWilman et al J 1 2003 . McGee 
et al., in prep), which are hard to reconcile with their model. On 
the other hand, the lower value of A/ trU nc~ 10 13 /!." 1 M that 
we advocate might have trouble e xplaining the large fr action of 
galaxies with [Oil] emission in the lPoggianti et all d2006h clusters 
at z ~ 0.5. Undoubtedly both models are greatly oversimplified 
and, moreover, there are important systematic uncertainties in the 
current data (especially in determining cluster masses and galaxy 
star formation rates) and statistical limitations resulting from small 
sample sizes. 

Similar constraints have also been deri ved in the pa s t from 
observations of radial gradients in clusters. iBalogh et al.l J2000l) 
used n-body simulations of the infall of substructure into clusters 
and concluded that, to match the radial gradients of star forma- 
tion rates, the star formation rates in cluster galaxies must decline 
on the timescale of a few Gyrs after entering the cluster. Signifi- 
cantly, they also found that the best match to radial gradients was 
provided if the star formation rate in the galaxy began to decline 
as soo n as it was found in a dark matter structure of group-size or 
larger. Ellingson et alj d200lh took this a step further and investi- 
gated the evolution of such gradients. They determined that 'field- 
like' galaxies b ecame early type galaxies on a 2-3 Gyr timescale. 
lEllingson et alj also inferred that if galaxies were transformed on 
the 3 Gyr timescale, than the galaxy infall rate into clusters be- 
tween z ~ 1.5 and z ~ 0.5 must have declined by ~ 20 %. Our 
results suggest that the infall rate of galaxies into clusters over the 
same span fell by ~ 15 %. This is a surprisingly good agreement 
given the large observational uncertainties at each step in this anal- 
ysis. 

Finally, it is instructive to reexamine the results of B09 in 
the context of our results. We have previously shown that B09 
finds a lower fraction of galaxies within groups and clusters than 
we do ( i|2.H . This is likely due to a subhalo completeness level 
which varies as a function of environment. They use the global 
number density of subhalos above their mass threshold and com- 
pare it against SDSS number densities to conclude that their global 
magnitude limit is ~ 0.3 L„ . However, when they compare num- 
ber densities of their subhalos within clusters with cluster obser- 
vations, they find that their cluster magnitude limit is ~ 0.5 L*. 
Using the red galaxy lumin osity function derived from a large sam- 
ple of galaxy clusters by iLu et alj d2009t) . a magnitude cut of 0.5 
L* instead of 0.3 L* reduces the number of cluster galaxies by ~ 
40 %. In other words, groups falling into their clusters could have 
~ 40 % fewer galaxies than would be expected from a consistent 
luminosity cut. Indeed, we find that this is on the order of the dis- 
crepancy between our results and those of B09. For instance, we 
have shown that, for 10 14 ' 2 hT 1 Mq clusters, ~ 35 % of galaxies 
have been accreted through 10 13 hT 1 Mq halos at z=0, while B09 



find only 24%. While this disagreement is significant for evaluat- 
ing the role of preprocessing in cluster assembly, a bigger factor is 
that the B09 clusters are not very massive. Indeed, these are smaller 
than the bulk of well studied clusters at intermediate and high red- 
shift. We have extended their analysis to more massive clusters and 
find, as B09 themselves anticipated, that group pre-processing is 
potentially much more important for more massive clusters. 

3.7 Towards a physically motivated model 

We have shown, by following the accretion of galaxies into groups 
and clusters, and making simple assumptions about the nature of 
environmental effects on galaxies, that the halo mass at which envi- 
ronmental effects begin to be induced on galaxies is approximately 
10 12 - 10 13 h' 1 Mq, and the time those effects take to manifest 
themselves is quite long (> 2 Gyr). Here, we address some of the 
more important simplifications we have made in constructing this 
model. 

The first simplification is that we have assumed that an en- 
vironmental effect will have a unique signature on the properties 
of galaxies. However, in comparing our model to, for instance, the 
fraction of red galaxies in clusters, we must acknowledge that there 
is more than one process which can make a galaxy red. In the local 
Universe, observations suggest that nearly all galaxies with stellar 
masses above 10 10 ft -1 Mq are red regardless of their environment 
jBaldrv et alj|2006h . However, as shown in Figure [10] in our sim- 
ple model the most massive cluster galaxies would still be red, a 
consequence of the fact that they have resided within massive dark 
haloes for a long time. This, combined with the fact that the more 
numerous low mass galaxies dominate the fraction of galaxies in a 
cluster, indicate that this is not a large complicating factor. 

Secondly, we have assumed that all galaxies display environ- 
mental effects after a s pecific time r tr u nc , reg ardless of their in- 
coming orbit. However. iMcCarfhv et alj d2008bh has shown in sim- 
ulations that the environmental effec t on an infalling galaxy is de- 
pendent on the orbit of that galaxy. iMcCarthv et alj also showed 
that the bulk of the environmental effect on an infalling galaxy 
occurs when the satellite is at its pericentre. The size of this ef- 
fect can be quantified by the variation in the time it takes a galaxy 
to fall from the virial radius to the pericentre of its orbit. In Fig- 
ure [TT] we show the distribution of times fo r a realistic di stribu- 
tion of infalling dark matter substructure from lBensonl J2005ri . ran- 
domly sampled 10,000 times. The distribution is shown as a mass- 
independent quantity, along with the best fit gaussian. A cluster of 
10 14 h,- 1 Mq has a R vir = 1.26 ft" 1 Mpc and V circ =400 km/s, 
which translates to a quite narrow distribution, with a dispersion of 
only ~ 0.2 Gyrs. This will not have significant implications for a 
timescale which is greater than 2 Gyrs. It is worth noting that the 
simulations of lBensor] were not adequate to quantify the effect of 
any host halo mass dependence of the orbital distribution, but the 
indications are that this will not have a significant impact for our 
purposes. Additionally, we have assumed that all galaxies entering 
a massive halo feel similar environmental effects, however, galax- 
ies with large pericentric distances may not feel strong effects, and 
thus predictions for the red fraction scatter and its mass dependence 
will still benefit from proper tracing of orbits in the future. 

Thirdly, we have assumed that the mass of a halo is the impor- 
tant quantity driving any environmental effects. In fact, most of the 
physical processes which could produce environmental effects are 
likely more sensitive to X-ray gas density or temperature (T x ). For 
galaxy clusters and massive groups the scatter in the M-T x relation 
is actually quite small (± 30 % at M = 10 14 ' 5 h' 1 Mq). However, 
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Figure 11. The distribution of times, t rp for an infalling dark matter sub- 
structure to reach its pericentre from the virial radius, R v i r , of a halo with 
circular velocity V c irc.vir- The bla c k line is the distribution of 10 5 ran- 
domly sampled orbits from iBensoiJ fe005l) and the dotted red line is the 
best fit gaussian. 



the Mass - X-ray Luminosity (M - L) relation, which is more sensi- 
tive to the gas density, does show significant sca tter at cluster mass 
scales I McC arthy et alj2004lBalogh et al.l2006h . But this scatter is 
driven by properties of the group and cluster cores, while at the ra- 
dius of a typical galaxy pericentre (0.2-0.3 times t he virial radius), 
the scatter from system to system is quite small ^McCarthy et al.l 
l2008ahlSun et al Jl 2008). So the bulk of ram pressure stripping will 
occur at radii where the gas density has little scatter from system to 
system. However, this analysis is limited to fairly massive groups 
and clsuters, as measurement of X-ray properties for typica l 10 12 
h^ 1 Mg haloes is quite difficult teeiprich & B6hringerl2002h . One 
theoretical indication of the size of this effect in low-mass groups 
comes from the scatter in the virial mass - circular vel ocity relation, 
which is approximately ± 15 % atM= 10 12 h" 1 Mm teullock et al.1 
1200 lh . The circular velocity is more indicative of the depth of the 
dark matter potential, and thus is likely more closely correlated 
with the gas density. Despite this, the size of this scatter is likely not 
a huge source of uncertainty in our model, given that we only make 
broad statements about the characteristic halo mass scale. Given all 
of these results, it is encouraging that our model does not appear too 
simple to give important insights to the behavior of environmental 
effects. 

The next step is to put this ad hoc model on a more physical 
basis. In particular, in our model we have specified that galaxies 
within a host halo are equally affected by environmental processes 
regardless of their position within the halo. But because the cooling 
rate of gas in a halo is density dependent, semi-analytic models treat 
galaxies at the center of halos (centrals) different from those not in 
the center (satellite s). While this distinction is still a simplification 
( Simha et al. 2008), we point out the similarity between a model 
with a A/trunc of 10 12 h^ 1 Mq and one where the environmental 
effect begins to occur when a galaxy becomes a satellite, as shown 
in Figure[8]. Our most favored Mtrunc model is essentially equiva- 
lent to choosing a physically motivated central/satellite model. 



We have also employed a fixed timescale for environmental ef- 
fects to occur. Ideally, we would like to link this timescale to a phys- 
ically motivated quantity, such as the orbital timescale of a galaxy 
in a cluster or group. This mean timescale is approximately con- 
stant for the groups and clusters in our mass range at a given red- 
shift epoch. However, because of the decreasing universe density 
with time, at high redshift the orbital timescale is actually smaller 
by a factor of ~ (l+z) 3,/2 . A timescale based on this would suggest 
that at z=1.5 the timescale is ~ 4 times shorter than the timescale at 
z=0. Unfortunately, directly implementing a timescale based on the 
orbital timescale would ignore several other complicating factors 
such as the evolution of cluster gas density profiles and the evolu- 
tion of galaxy sizes and densities. Exploring these issues in a full 
semi-analytic galaxy model is the important next step forward. 



4 DISCUSSION AND CONCLUSIONS 

We have used the stellar mass and merger trees produced by the 
semi-analytic galaxy catalogues of F08 to follow the accretion of 
galaxies into groups and clusters at four different redshift epochs 
(z=0,0.5,1.0 and 1.5) for samples of galaxies with stellar mass 
M > 10 9 /i _1 Mq. By tracking galaxies through the hierarchy of 
structure formation we are able to examine the effect that environ- 
mental processes may have on the galaxy population of groups and 
clusters. Further, by adopting a simple model for the environmental 
effects, we are able to make strong claims about the timescale and 
mass threshold on which environmental effects occur. Our main re- 
sults are summarized as follows: 

• Clusters at all redshifts examined have had a significant frac- 
tion of their galaxies accreted through galaxy groups. For instance, 
10 14 ' 5 h' 1 Mq mass clusters at z=0 have had ~ 40% of their 
galaxies (M stc n ar > 10 9 h -1 Mq) accreted through halos with 
masses greater than 10 13 h~ x Mq. At higher redshifts fewer galax- 
ies are accreted through massive halos. Only ~ 25 % of galaxies 
have been accreted through 10 13 h' 1 M Q into 10 14 5 h" 1 M Q 
mass clusters at z=1.5. 

• We find only a moderate difference in the stellar mass accre- 
tion history and the galaxy accretion history at high cluster mass. 
That is, more massive galaxies are accreted preferentially through 
groups. While 45% of galaxies in 10 15 h' 1 M mass clusters 
at z=0 are accreted through halos with masses greater than 10 13 
h~ Mq, 50% of the stellar mass is ac creted through the sa me halo 
mass range. Contrary to the study of iBerrier et alj ( 120091) . we do 
not see a large difference between the galaxy assembly of clusters 
and the mass assembly of clusters. 

• Following from the previous point, we find that the extent to 
which galaxies are pre-processed in groups before falling into clus- 
ters depends on the stellar mass of the infalling galaxy. For a 10 14 ' 5 
ir 1 Mq mass cluster, 73 % of galaxies with stellar masses greater 
than 10 10 ' 5 ft -1 M are accreted through 10 12 h^ 1 Mq systems, 
while only 50% of 10 9 to 10 10 h' 1 Mq are accreted through the 
same systems. Further, we find that in the accretion through group 
sized halos increases at late times when compared to the accretion 
into the cluster during early times. 

• We have shown that the fraction of isolated galaxies infalling 
into z=0 groups and clusters is remarkably independent of the final 
cluster mass. 5-6 % of the final cluster galaxies are accreted per 
Gyr for the last 10 Gyrs. Thus if a galaxy begins to be affected by 
its environment soon after becoming a satellite galaxy, and the time 
it takes for that effect to manifest itself is constant with halo mass, 



The accretion of galaxies into groups and clusters 13 



then a similar fraction of galaxies are affected in each cluster above 
a halo mass of 10 13 ft" 1 M Q . 

• Despite the previous result, observing a cluster of the same 
halo mass at each redshift epoch implies different accretion rates 
of isolated galaxies, from 5-6 % per Gyr at z=0 to 15% per Gyr at 
z=1.5. Thus, in effect, the Butcher Oemler effect may be qualita- 
tively explained by the shorter time available for cluster assembly 
at higher redshift. 

• We find that combining the simple observations of the exis- 
tence of a significant Butcher Oemler effect at z=0.5 and the obser- 
vations that galaxies within groups display significant environmen- 
tal effects with galaxy accretion histories justifies striking conclu- 
sions. Namely, that the dominant environmental process must begin 
to occur in halos of 10 12 - 10 13 h -1 Mq and act over timescales 
of > 2 Gyrs. This supports a long lifetime, gentle mechanism like 
strangulation. 

• This simple model predicts that by z=1.5 galaxy groups and 
clusters will display little to no environmental effects. This con- 
clusion may have limit the effectiveness of red sequence cluster 
finding methods at high redshift. 

In essence, we have seen that systematic observations of in- 
termediate and high redshift clusters and groups have the power 
to strongly constrain the mechanisms which induce environmental 
transformations on galaxies. However, because of the significant 
cluster to cluster variations in environmental effects, it is impor- 
tant that the method for selecting galaxy clusters and groups for 
observation must be easily and accurately reproducible in cosmo- 
logical simulations. Only this will allow the careful testing of mod- 
els against observations. In a future paper we will compare the best 
available data on groups and clusters at a variety of redshift epochs 
to further constrain the dominant environmental processes. 

Significant progress on the implications of strangulation and 
the physical processes involved will need more extensive hydrody- 
namical simulations. The simulations of ram pressure stripp ing of 
the hot haloes of infalling galaxies by McCarthy et al. (20083) is a 
significant step forward. However, there are important unknowns. 
In particular, how effective are low mass group halos in stripping 
the infalling galaxies? Unfortunately, this is sensitively dependent 
on how the gas is distributed in both the infalling galaxy and the 
group halos. Indeed, the effectiveness of strangulation is also de- 
pendent on the strength of star formation feedback, and how re- 
heated galaxy gas is distributed and stripped from the galaxy. The 
behavior of galaxies within small groups which subsequently fall 
into massive clusters is also unclear. To what extent are galaxies 
"shielded" by their local group from further gas stripping? Encour- 
agingly, large scale hydrodynamical simulatio ns are beginning to 
be able to address some of these questions (e.g. lCrain et alj |2009). 

So, while there is much room for improvement in understand- 
ing the details of galaxy - environment interactions, our results 
have shown that the galaxy accretion histories of groups and clus- 
ters combined with a simple model strongly suggest that the domi- 
nant environmental effect occurs over long time scales and is effec- 
tive in low mass halos. In a future paper, we will examine these in- 
sights by making a quantitative comparison between semi-analytic 
models and the best available cluster, group and field data to z ~ 1 . 
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APPENDIX: FULL ACCRETION HISTORIES 

Here we show the complete accretion histories for each bin of clus- 
ter mass and for all four redshift epochs. They are presented for 
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both the galaxy accretion (Figure [72t and for stellar mass accre- 
tion (Figure[T3t. Figure[T2lshows the cumulative distribution of ac- 
creted cluster galaxies which reside in a host halo mass of a given 
size prior to accretion into the final cluster at each of four epochs of 
observation. Because galaxies are on average more massive in more 
massive halos, this accretion history does not agree completely with 
dark matter accretion histories. Therefore, we present the complete 
stellar mass accretion histories in Figure[T3] Again, this shows the 
cumulative distribution of the accreted stellar mass as a function of 
the galaxy's host halo mass at the time of accretion. 
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Figure 12. The cumulative distribution of cluster galaxies which reside in a host halo of a given mass at the time of accretion into the final cluster halo. In the 
left panel is the accretion history of 8 composite clusters of a given final host mass at z=0, while in the left middle (right middle) [right] panel is a separate 
final cluster sample at z=0.5 (z=l) [z=1.5]. All cluster galaxies have final stellar masses of M > 10 s Mq. The mass range bins were defined in Table 1, and 
are shown for all bins containing more than one cluster. 
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Figure 13. The cumulative distribution of accreted stellar mass which reside in a host halo of a given size at the time of accretion into the final cluster halo. In 
the left panel, is the accretion history of 8 composite clusters of a given final host mass at z=0. The left middle (right middle ) [right] panel is for a separate 
final cluster sample at z=0.5 (z=l) [z=1.5]. All cluster galaxies have final stellar masses of M > 1O 9 M . The mass range bins were defined in Table 1 , and 
are shown for all bins containing more than one cluster. 



